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1 General remarks 

Any object with a net positive charge on one end and a net negative charge 
on the other end possesses an electric dipole moment. In ordinary classical 
electromagnetism this dipole moment is a vector quantity that can point in any 
direction, and is subject to electrical forces that are fairly straightforward to 
formulate mathematically. However, for a quantum mechanical object like an 
atom or molecule, the strength and orientation of the object's dipole moment 
can depend strongly on the object's quantum mechanical state. This is a subject 
that becomes relevant in low temperature molecular samples, where an ensemble 
of molecules can be prepared in a single internal state, as described elsewhere 
in this volume. In such a case, the mathematical description becomes more 
elaborate, and indeed the dipole-dipole interaction need not take the classical 
form given in textbooks. The description of this interaction is the subject of 
this chapter. 

We approach this task in three steps: first, we introduce the ideas of how 
dipoles arise in quantum mechanical objects; second, we present a formalism 
within which to describe these dipoles; and third, we give examples of the for- 
malism that illustrate some of the basic physics that emerges. The discussion 
will explore the possible energy states of the dipole, the field generated by the 
dipole, and the interaction of the dipole with another dipole. We restrict the 
discussion to a particular "minimal realistic model," so that the most important 
physics is incorporated, but the arithmetic is not overwhelming. 

Although we discuss molecular dipoles in several contexts, our main focus 
is on polar molecules that possess a A doublet in their ground state. These 
molecules are the most likely, among diatomic molecules at least, to exhibit their 
dipolar character at moderate laboratory field strengths. A-doubled molecules 
have another peculiar feature, namely, their ground states possess a degeneracy 
even in an electric field. This means that there is more than one way for such 
a molecule to align with the field; the two possibilities are characterized by 
different angular momentum quantum numbers. This degeneracy leads to novel 
properties of both the orientation of a single molecule's dipole moment and the 
interaction between dipoles. In the examples we present, we focus on revealing 
these novel features. 



We assume the reader has a good background in undergraduate quantum me- 
chanics and electrostatics. In particular, the ideas of matrix mechanics, Dirac 
notation, and time-independent perturbation theory are used frequently. In ad- 
dition, the reader should have a passing familiarity with electric dipoles and 
their interactions with fields and with each other. Finally, we will draw heav- 
ily on the mathematical theory of angular momentum as applied to quantum 
mechanics, as described in Appendix A of this volume, and in more detail in 
the classic treatise of Brink and Satchler [I]. When necessary, details of the 
structure of diatomic molecules have been drawn from Brown and Carrington's 
recent authoritative text [5]. 

2 Review of Classical Dipoles 

The behavior of a polar molecule is largely determined by its response to electric 
fields. Classically, an electric dipole appears when a molecule has a little bit of 
positive charge displaced a distance from a little bit of negative charge. The 
dipole moment is then a vector quantity that characterizes the direction and 
magnitude of this displacement: 

where the charge is displaced from a particular origin. Because we 
are interested in forces exerted on molecules, we will take this origin to be the 
center of mass of the molecule. (Defining fl — Q would instead identify the 
center of charge of the molecule - quite a different thing!) By convention, the 
dipole moment vector points away from the negative charges, and toward the 
positive charges, inside the molecule. 

A molecule has many charges in it, and they are distributed in a complex way, 
as governed by the quantum mechanical state of the molecule. In general there 
is much more information about the electrostatic properties of the molecule 
than is contained in its dipole moment. However, at distances far from the 
molecule (as compared to the molecule's size), these details do not matter. The 
forces that one molecule exerts on another in this limit is strongly dominated 
by the dipole moments of the two molecules. We consider in this chapter only 
electrically neutral molecules, so that the Coulomb force between molecules is 
absent. In this limit, too, the details of the dipole moment's origin are irrelevant, 
and we consider the molecule to be a "point dipole," whose dipole moment is 
characterized by a magnitude /i and a direction /t. 

If a dipole fl is immersed in an electric field £, its energy depends on the 
relative orientation of the field and the dipole, via 

Eel = -ft - £■ 

This follows simply from the fact that the positive charges will be pulled in the 
direction of the field, while the negative charges are pulled the other way. Thus 
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a dipole pointing in the same direction as the field is lower in energy than a 
dipole pointing in exactly the opposite direction. In classical electrostatics, the 
energy can continuously vary between these two extreme limits. 

As an object containing charge, a dipole generates an electric field, which is 
given, as usual, by the gradient of an electrostatic potential, ^molecule = — V$(r). 
For a point dipole the potential $ is given by 

m = (1) 

where f = rf denotes the point in space, relative to the dipole, at which the field 
is to be evaluated ^Sj • The dot product in ^ gives the field <& a strong angular 
dependence. For this reason, it is convenient to use spherical coordinates to 
describe the physics of dipoles, since they explicitly record directions. Setting 
f = {6, (p) and /t — (a, [3) in spherical coordinates, the dipole potential becomes 

<& = (cos a cos + sinasin0cos(/3 — (j))) . (2) 

For the most familiar case of a dipole aligned along the positive z-axis (a = 0), 
this yields the familiar result $ = ji cos 6 /r^. This potential is maximal along 
the dipole's axis (0 = or tt), and vanishes in the direction perpendicular to 
the dipole {0 = 7r/2). 

From the results above we can evaluate the interaction potential between 
two dipoles. One of the dipoles generates an electric field, which acts on the 
other. Taking the scalar product of one dipole moment with the gradient of the 
dipole potential ^ due to the other, we obtain 3J 

T/ ^ m Ml • A*2 - 3(^1 • R){fi2 ■ R) 
Vd{R) ^3 , 

where R — RR is the relative coordinate of the dipoles. This result is general 
for any orientation of each dipole, and for any relative position of the pair of 
dipoles. In a special case where both dipoles ar aligned along the positive z-axis, 
and where the vector connecting the centers-of-mass of the two dipoles makes 
an angle with this axis, the dipole-dipole interaction takes a simpler form: 

T . / ^^ 1 — 3 cos^ 9 

Vd{r) ^^Xl^i2 . (3) 

Note that the angle 9 as used here has a different meaning from the one in Eq. 

We will use 9 is both contexts throughout this chapter, hopefully without 
causing undue confusion. The form ^ of the interaction is useful for illustrating 
the most basic fact of the dipolc-dipolc interaction: if the two dipoles line up 
in a head-to-tail orientation [9 — or tt), then < and they attract one 
another; whereas if they lie side- by-side [9 = tt/2), then > and they repel 
one another^ 

^This expression ignores a contact potential that must be associated to a point dipole 
to conserve lines of electric flux 3 . However, real molecules are not point dipoles, and the 
electrostatic potential differs greatly from this dipolar form at length scales inside the molecule, 
scales that do not concern us here. 
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Our main goal in this chapter is to investigate how these classical results 
change when the dipoles belong to molecules that are governed by quantum 
mechanics. In Sec. 3 we evaluate the energy of a dipole exposed to an external 
field; in Sec. 4 we consider the field produced by a quantum mechanical dipole; 
and in Sec. 5 we address the interaction between two dipolar molecules. 

3 Quantum mechanical dipoles in fields 

Whereas the classical energy of a dipole in a field can take a continuum of values 
between its minimum and maximum, this is no longer the case for a quantum 
mechanical molecule. In this section we will establish the spectrum of a polar 
molecule in an electric field, building from a set of simple examples. To start, we 
will define the laboratory z-axis to coincide with the direction of an externally 
applied electric field, so that 8 = Sz. In this case, the projection of the total 
angular momentum on the z-axis is a conserved quantity. 

3.1 Atoms 

Our main focus in this chapter will be on electrically polarizable dipolar molecules. 
But before discussing this in detail, we first consider the simpler case of an 
electrically polarizable atom, namely, hydrogen. This will introduce both the 
basic physics ideas, and the angular momentum techniques that we will use. In 
this case a negatively charged electron separated a distance r from a positively 
charged proton forms a dipole moment jl = —ef. 

Because dipoles require us to consider directions, it is useful to cast the unit 
vector f into its spherical components [1]: 

(-±z^) =T%/2Ci±i(0,0), i=C7io(0,0). 
\r r / r 

Here the C"s are reduced spherical harmonics, related to the usual spherical 
harmonics by [T] 

and given explicitly for fc = I by (Appendix A) 

Ci±i{e(t>) = T^sinee^'^ Cio(6l(/>) = cos 61. (4) 
v2 

In general, it is convenient to represent interaction potentials in terms of the 
functions Ckq (since they do not carry extra factors of 47r), and to use the 
functions Y^q as wave functions in angular degrees of freedom (since they are 
already properly normalized, by (Y/m|y/'m') = ^ii'^mm')- Integrals involving 
the reduced spherical harmonics are conveniently related to the 3-j symbols of 
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angular momentum theory, for example: 



J d(cos e)dcpck,g, mck,,, mck.q, m 

^ A f ^2 ks \ I' ki k2 ks \ 
~ qi q2 ks J \ )■ 

The 3-j symbols, in parentheses, are related to the Clebsch-Gordan coefficients. 
They are widely tabulated and easily computed for applications (Appendix A). 
In terms of these functions, the Hamiltonian for the atom-field interaction is 

iJci = -i-ef) ■ £ = ez£ = ercos(6l)f = er£Cw{e). (5) 

The possible energies for a dipole in a field are given by the eigenvalues of the 
Hamiltonian ([5]). To evaluate these energies in quantum mechanics, we identify 
the usual basis set of hydrogenic wave functions, \nlm), where we ignore spin 
for this simple illustration: 

{r,e,4>\nlm) ^ fni{r)Yi„,{e,cp). 

The matrix elements between any two hydrogenic states are 

{nlm\- p- £\n'l'm') ^ {er)£ J d{cose)d(l)Yi*^CioYi>„i' (6) 

-ier)£vwTWFTT)^^^ I ;;)( I j;) 

Here (er) = J r'^drfni{r)rfn'i' is an effective magnitude of the dipole moment, 
which can be analytically evaluated for hydrogen [4J. 

Some important physics is embodied in Eqn. ([6]). First, the electric field 
defines an axis of rotational symmetry (here the z axis). On general grounds, 
we therefore expect that the projection of the total angular momentum of the 
molecule onto this axis is a constant. And indeed, this is built into the 3-j 
symbols: since the sum of all m quantum numbers in a 3-j symbol should add 
to zero, Eqn. ([6]) asserts that m — m' , and the electric field cannot couple two 
different to's together. 

A second feature embodied in ^ is the action of parity. The hydrogenic 
wavefunctions have a definite parity, i.e., they either change sign, or else remain 
invariant, upon converting from a coordinate system {x,y,z) to a coordinate 
system {—x, —y, —z). The sign of the parity-changed wave function is given by 
(—1)'. Thus an s-state (/ = 0) has even parity, while a p-state {I — 1) has odd 
parity. For an electric field pointing in a particular direction, the Hamiltonian 
([5]) itself has odd parity, and thus serves to change the parity of the atom. For 
example, it can couple the s and p states to each other, but not to themselves. 
This is expressed in ^ by the first 3-j symbol, whose symmetry properties 
require that / + I + / = even. This seemingly innocuous statement is the 
fundamental fact of electric dipole moments of atoms and molecules. It says 
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that, for example the Is ground state of hydrogen, with I = I' = 0, does not, by 
itself, respond to an electric field at all. Rather, it requires an admixture of a p 
state to develop a dipole moment 

To evaluate the influence of an electric field on hydrogen, therefore, we must 
consider at least the nearest state of opposite parity, which is the 2p state. These 
two states are separated in energy by an amount Eis2p- Considering only these 
two states, and ignoring any spin structure, the atom-plus-field Hamiltonian is 
represented by a simple 2x2 matrix: 



where the dipole matrix element is given by the convenient shorthand /x = 
(ls,m = 0|ez|2p,m = 0) = 128%/2eao/243 g]. Of course there are many more 
p states that the Is state is coupled to. Plus, all states are further complicated 
by the spin of the electron and (in hydrogen) the nucleus. Matrix elements for 
all these can be constructed, and the full matrix diagonalized to approximate 
the energies to any desired degree of accuracy. However, we are interested here 
in the qualitative features of dipoles, and so limit ourselves to Eqf.(l7|). 
The Stark energies are thus given approximately by 



This expression illustrates the basic physics of the quantum mechanical dipole. 
First, there are necessarily two states (or more) involved. One state decreases 
in energy as the field is turned on, representing the "normal" case where the 
electron moves to negative z and the electric dipole moment aligns with the field. 
The other state, however, increases in energy with increasing field and represents 
the dipole moment anti-aligning with the field. Classically it is of course possible 
to align the dipole against the field in a state of unstable equilibrium. Similarly, 
in quantum mechanics this is a legitimate energy eigenstate, and the dipole will 
remain anti-aligned with the field in the absence of perturbations. 

A second observation about the energies ([5]) is that the energy is a quadratic 
function of E at low field, and only becomes linearly proportional to £ at higher 
fields. Thus the permanent dipole moment of the atom, defined by the zero-field 
limit 



vanishes. The atom, in an energy eigenstate in zero field, has no permanent 
electric dipole moment. This makes sense since, in zero field, the electron's 
position is randomly varied about the atom, lying as much on one side of the 
nucleus as on the opposite side. 

■^These remarks are not strictly true. The ground state of hydrogen aheady has a small 
admixture of odd-parity states, due to the parity-violating part of the electroweak force. This 
effect is far too small to concern us here, however. 






(8) 



/^permanent — 
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The transition from quadratic to linear Stark effect is an example of a com- 
petition between two tendencies. At low field, the dominant energy scale is the 
energy splitting Eis2p between opposite parity states. At higher field values, 
the interaction energy with the electric field becomes stronger, and the dipole 
is aligned. The value of the field where this transition occurs is found roughly 
by setting these energies equal to find a "critical" electric field: 

^critical = -£'ls2p/2(i. 

For atomic hydrogen, this field is on the order of 10^ V/cm. However, at this 
field it is already a bad approximation to ignore that fact that there are both 
2pi/2 and 2p3/2 states, as well as higher- lying p states, and further coupling 
between p, d, etc., states. We will not pursue this subject further here. 

3.2 Rotating molecules 

With these basics in mind, we can move on to molecules. We focus here on di- 
atomic, heteronuclear molecules, although the principles are more general. We 
will consider only electric fields so small that the electrons cannot be polarized 
in the sense of the previous section; thus we consider only a single electronic 
state. However, the charge separation between the two atoms produces an elec- 
tric dipole moment jl in the rotating frame of the molecule. We assume that 
the molecule is a rigid rotor and we will not consider explicitly the vibrational 
motion of the molecule, focusing instead solely on the molecular rotation. (More 
precisely, we consider jJ, to incorporate an averaging over the vibrational coordi- 
nate of the molecule, much as the electron-proton distance r was averaged over 
for the hydrogen atom in the previous section.) 

As a mathematical preliminary, we note the following. To deal with molecules, 
we are required to transform freely between the laboratory reference frame and 
the body-fixed frame that rotates with the molecule. The rotation from the lab 
frame (x, y, z) to the body frame (x', y', z') is governed by a set of Euler angles 
(a, /3, 7) (Appendix A). The first two angles a = 4>, (3 = 9 coincide with the 
spherical coordinates {9, (j)) of the body frame's z' axis. By convention, we take 
the positive z' direction to be parallel to the dipole moment pi. The third Euler 
angle 7 serves to orient the x' axis in a desired orientation within the body 
frame; it is thus the azimuthal angle of rotation about the molecular axis. 

Consider a given angular momentum state |jm) referred to the lab frame. 
This state is only a state of good m in the lab frame, in general. In the body 
frame, which points in some other direction, the same state will be a linear 
superposition of different m's, which we denote in the body frame as cij's to 
distinguish them. Moreover, this linear superposition will be a function of the 
Euler angles, with a transformation that is conventionally denoted by the letter 
D: 

D{aP^)\juj) = ^ \ jm){jm\D{a, P,^)\juj) 

m 
m 
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This last line defines the Wigner rotation matrices, whose properties are widely 
tabulated. For each j, Di^^ is a unitary transformation matrix; note that 
a rotation can only change m-type quantum numbers, not the total angular 
momentum j . One of the more useful properties of the D matrices, for us, is 

dadcos(/3)d7 Di^,^MPl)Di\^.^o.Pl)D'^,.,MPl) = 

( Ji n J3 \ ( n n n \ 



TOi 7712 "^-3 / V ^1 ^2 777(^3 

Because the dipole is aligned along the molecular axis, and because the 
molecular axis is tilted at an angle (i with respect to the field, and because 
the field defines the z-axis, the dipole moment is defined by its magnitude [i 
times a unit vector with polar coordinates a) . The Hamiltonian for the 
molecule-field interaction is given by 

i/ei = ■ £ = -ti£Cioil3a) = -fxSDi^liap-f). (10) 

For use below, we have taken the liberty of rewriting Cio as a _D-function; since 
the second index of D is zero, this function does not actually depend on 7, so 
introducing this variable is not as drastic as it seems. 

To evaluate energies in quantum mechanics we need to choose a basis set 
and take matrix elements. The Wigner rotation matrices are the quantum 
mechanical eigenfunctions of the rigid rotor. With normalization, these wave 
functions are 



(a/?7|7i77i„A„) = y^^^^":„A„ 



As we did for hydrogen, we here ignore spin. Thus 77 is the quantum number of 
rotation of the atoms about their center of mass, 777„ is the projection of this 
angular momentum in the lab frame, and A„ is its projection in the body frame. 
In this basis, the matrix elements of the Stark interaction are computed using 
dH), to yield 

(n77i„A„| - fici ■ £\n 777„A ) (11) 



-Alelf (-l)™"-^'y(271+ 1)(277' + 1) 



77 l77\/77 I77 

-777„ m'^ I \ -A„ A^j 



In an important special case, the molecule is in a S state, meaning that the 
electronic angular momentum projection A„ = 0. In this case, Eqn. (jlip reduces 
to the same expression as that for hydrogen, apart from the radial integral. This 
is as it should be: in both objects, there is simply a positive charge at one end 
and a negative charge at the other. It does not matter if one of these is an 
electron, rather than an atom. More generally, however, when An ^ there will 
be a complicating effect of lambda-doubling, which we will discuss in the next 
section. 
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Thus the physics of the rotating dipole is much the same as that of the 
hydrogen atom. Eqn. pT|) also asserts that, for a S state with A„ = 0, the 
electric field interaction vanishes unless n and n' have opposite parity. For such 
a state, the parity is related to the parity of n itself. Thus, for the ground state 
of a molecule with n = 0, the electric field only has an effect by mixing this 
state with the next rotational state with n' = 1. These states are split by an 
energy i^iot = 2i?e, where Be is the rotational constant of the moleculeU 

We can formulate a simple 2x2 matrix describing this situation, as we did 
for hydrogen: 



where the dipole matrix element is given by the convenient shorthand notation 
^ ~ {nmnO\fiq^o\n'mnO) . There is one such matrix for each value of m„. Of 
course there are many more rotational states that these states are coupled to. 
Plus, all states would further be complicated by the spins (if any) of the electrons 
and nuclei. 

The matrix (fT2|) can be diagonalized just as ([7|) was above, and the same 
physical conclusions apply. Namely, the molecule in a given rotational state 
has no permanent electric dipole moment, even though there is a separation of 
charges in the body frame of the molecule. Second, the Stark effect is quadratic 
for low fields, and linear only at higher fields, with the transition occurring at 
a "critical field" 



To take an example, the NH molecule posesses a '^E ground state. For this state, 
ignoring spin, the critical field is of the order 7x10^ V/cm. This is far smaller 
than the field required to polarize electrons in an atom or molecule, but still 
large for laboratory-strength electric fields. Diatomic molecules with smaller 
rotational constants, such as LiF, would have correspondingly smaller critical 
fields. In any event, by the time the critical field is applied, it is already a bad 
approximation to ignore coupling to the other rotational states of the molecule, 
which must be included for an accurate treatment. We do not consider this 
topic further here. 

3.3 Molecules with lambda-doubling 

As we have made clear in the previous two sections, the effect of an electric 
field on a quantum mechanical object is to couple states of opposite parity. For 
a molecule in a 11 or A state, there are often two such parity states that are 
much closer together in energy than the rotational spacing. The two states are 
said to be the components of a "A-doublet." Because they are close together 
in energy, these two states can then be mixed at much smaller fields than are 
required to mix rotational levels. The physics underlying the lambda doublet is 
rather complex, and we refer the reader to the literature for details [2 [5] . 

■^In zero field, the state with rotational quantum number n has energy Be.n{n + 1). 




(12) 



i?rot/2/i. 



(13) 
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However, in broad terms, the argument is something hke this: a 11 state has 
an electronic angular momentum projection of magnitude 1 about the molecular 
axis. This angular momentum comes in two projections, for the two sense of ro- 
tation about the axis, and these projections are nominally degenerate in energy. 
The rotation of the molecule, however, can break the degeneracy between these 
levels, and (it so happens) the resulting nondegenerate eigenfunctions are also 
eigenfunctions of parity. The main point is that the resulting energy splitting is 
usually quite small, and these parity states can be mixed in fields much smaller 
than those required to mix rotational states. 

To this end, we modify the rigid-rotor wave function of the molecule to 
incorporate the electronic angular momentum: 



{a(3j\jm^) = ^I^^D^^^{apj). (14) 



Here j is the total (rotation-plus-electronic) angular momentum of the molecule, 
and m and w are the projections of j on the laboratory and body-fixed axes, 
respectively. Using the total j angular momentum, rather than just the molec- 
ular rotation n, marks the use of a "Hund's case a" representation, rather than 
the Hund's case b that was implicit in the previous section (Ref. [5]; see also 
Appendix B of this volume). 

In this basis the matrix element of the electric field Hamiltonian (fTOl) becomes 



(jmtj| - flci ■ £\j'm'uj') (15) 
= -Aicif(-l)'"-'^(2j + l) ' ^ ^ f \ l J 1 / 



-m m J \ —ui Lo 

In (fT5|l . the 3-j symbols denote conservation laws. The first asserts that m = m! 
is conserved, as we already knew. The second i-j symbol adds to this the fact 
that u) — Ld' . This is the statement that the electric field cannot exert a torque 
around the axis of the dipole moment itself. Moreover, in the present model we 
assert that j = j', since the next higher-lying j level is far away in energy, and 
only weakly mixes with the ground state j. With these approximations, the 
3-j symbols have simple algebraic expressions, and we can simplify the matrix 
element: 

{jmuo\ — /id • £\imuj) — -—fiS- 



The physical content of this expression is illustrated in Fig. 1. Notice that both 
m and lu can have a sign, and that whether the energy is positive or negative 
depends on both signs. 

An essential point is that there is not a unique state representing the dipole 
aligned with the field. Rather, there are two such states, distinguished by 
different angular momentum quantum numbers but possessing the same energy. 
To distinguish these in the following we will refer to the two states in figures 
1(a) and 1(b) as molecules of type \a) and type \b), respectively. Likewise, for 
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molecules nominally anti- aligned with the field, we will refer to types |c) and 
Id), corresponding to the two states in figures 1(c) and 1(d). The existence of 
these degeneracies will lead to novel phenomena in these kinds of molecules, as 
we discuss below. 

As to the lambda doubling, it is, as we have asserted, diagonal in a basis 
where parity is a good quantum number. In terms of the basis wave 
functions of well-defined parity are given by the linear combinations 

\jmuje) = [\jmoj) + e\jTnui)] . (16) 

Here we define Cu = \u)\, the absolute value of lu. For a given value of m, 
the linear combinations of ±(D in p6p are distinguished by the parity quantum 
number e = ±1. It is straightforward to show that in the parity basis (fTB|) the 
A-doubling is off-diagonal. 

The net result is that for each value of m, the Hamiltonian for our lambda- 
doubled molecule can be represented as a two-by-two matrix, similar to the ones 
above: 



A/2 Q 

where A is the lambda doubling energy, i.e., the energy difference between the 
two parity states, and 

Q = fit- 



is a manifestly positive quantity. This is the Hamiltonian we will treat in the 
remainder of this chapter. The difference here from the previous subsections is 
that the basis is now , which diagonalizes the electric field interaction rather 
than the zero-field Hamiltonian. This change reflects our emphasis on molecules 
in strong fields where their dipole moments are made manifest. The zero- field 
A-doubling Hamiltonian is considered, for the most part, to be a perturbation. 

The mixing of the strong field states due to the A-doubling interaction is 
conveniently given by a mixing angle S„n which we define as follows: 

\muje = +) = COS Sm\mLo) + sin Smlm — uj) 

\muje ^ —) = — sin Smlmuj) + cos Smlm — Co) (17) 

Explicitly, the mixing angle as a function of field is given by 

A/2 

" n I ^ /r-r^ - -tan<5_|™|, 

in terms of the energy Q and the dimensionless parameter 

A 

2Q' 
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Figure 1: Energetics of a polar molecule in an electric field. The molecule's 
dipole moment points from the negatively charged atom (large circle) to the 
positively charged atom (small circle) as indicated by the thick arrow. The 
dashed line indicates the positive direction of the molecules body axis, while 
the vertical arrow represents the direction of the applied electric field. The 
dipole aligns with the field, on average, if either i) the angular momentum j 
aligns with the field and m > 0, cj > (a); or ii) j aligns against the field and 
m < 0, a; < (b). Similar remarks apply to dipoles that anti- align with the 
field (c,d). 
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Notice that with this definition, Sm is positive when m is positive, and 5^„i = 
—6m- The energies of these states are conveniently summarized by the expres- 
sion 

JU + 1) 

This is a very compact way of writing the resuhs that will facilitate writing the 
expressions below. Notice that the intuition afforded by Figure 1 is still intact in 
these energies, but by replacing the sign of e for the sign oiuj. Thus states with 
me > have negative energy, while states with me < have positive energy. 

The case where m = must be handled slightly differently, since in this case 
the electric field energies Q — 0. We can still write the eigenstates in the form 
p7|) . provided that we set 

and understand that the corresponding energies, independent of field, are simply 

Eo^e - -eA/2. 

Like the dipoles considered above, this model has a quadratic Stark effect 
at low energies, rolling over to linear at electric fields exceeding the critical field 
given by setting Q — A/2. This criterion gives a critical field of 

_ Aj(j + 1) , . 

To take an example, consider the ground state of OH, which has j ~ lu ~ 3/2, 
fi = 1.7 Db, and a A-doublet splitting of 0.06 cm~^. In its m = 3/2 ground 
state, its critical field is ~ 1600 V/cm. Again, we have explicitly ignored the spin 
of the electron. The parity states are therefore easily mixed in fields that are 
both small enough to easily obtain in the laboratory, and small enough that no 
second-order coupling to rotational or electronic states needs to be considered. 
Keeping the relatively small number of molecular states is therefore a reasonable 
approximation, and highly desirable as it simplifies our discussion|f| 



4 The field due to a dipole 

Once polarized, each of the A-doubled molecules discussed above is itself the 
source of an electric field. The field due to a molecule is given above by Eqn. ([T]). 
In what follows, it is convenient to cast this potential in terms of the spherical 
tensors defined above and in Appendix A: 

$ (r) = ^ = ^ ^(-l)«Ci,(a/3)Ci_,(0,/>), (19) 

9 

OH there is also the u> = 1/2 state to consider, but it is also far away in energy as 
compared to the A-doublet, and we ignore it. It does play a role in the fine structure of OH, 
however, and this should be included in a quantitative model of OH. 
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That this form is correct can be verified by simple substitution, using the defi- 
nitions (HI . It seems at first unnecessarily complicated to write (|19p in this way. 
However, the effort required to do so will be rewarded when we need to evaluate 
the potential for quantum mechanical dipoles below. 

For a classical dipole, defining the direction (a/3) of its dipole moment would 
immediately specify the electrostatic potential it generates according to Eq. 
(|19p . However, in quantum mechanics the potnetial will result from suitably 
averaging the orientation of /j, over the distribution of {a/]) dictated by the 
molecule's wave function. To evaluate this, we need to evaluate matrix elements 
of (fTO|) in the basis of energy eigenstates P?]) . This is easily done using Eq.®, 
along with formulas that simplify the 3-j symbols. The result is 



{jmuj\Cig\jm'Lj') =5^^, " { rn, g = S, (20) 



j(i + l) 



{j-rn){j+m' ) 

2 ' y 



where the quantum numbers m and lo are signed quantities. Also note that 
this integration is over the molecular degrees of freedom (a/3) in lfT9|) . This still 
leaves the angular dependence on (0(/>), which characterizes the field in space 
around the dipole. 

From expression (PO)) it is clear that this matrix element changes sign upon 
either i) reversing the sign of both m and m! or ii) changing the sign of ui. 
Moreover, {jmuj\Ciq\j'm' — lo) =0, since to is conserved by the electric field. 
Using these observations, we can readily compute the matrix elements of $ in 
the dressed basis (|16p . Generally they take the form 



(jmwe|$(f)|jmW) = (jm(:;e|Ci,|jmW)(-l)*^Ci_,(00). (21) 

The matrix elements in front of (j2ip represent the quantum mechanical man- 
ifestation of the dipole's orientation. These matrix elements follow from the 
above definition of the dressed states, (|17p . Explicitly, 

(jmu),e|Cig|jm'w, e) = ecoii(5ra + 5m'){imQ\Ciq\jm' u)) 
{jmu},-\Ciq\jm'u},+) = {jmu),+\Ciq\jm'u),-') (22) 
= - sin((5„i + 5,n'){3mu}\Ciq\jm' u)) . 

In all these expressions, the value of q is set by angular momentum conservation 
to q = m — m'. 

This description, while complete within our model, nevertheless remains 
somewhat opaque. Let us therefore specialize it to the case of a particular 
energy eigenstate \jmuje). In this state, the averaged electrostatic potential of 
the dipole is 

f meui „^ \ cos 9 . . 

(*(^)) = ('^7(7TT) ^^^^ 
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Here the factor in parentheses is a quantum mechanical correction to the mag- 
nitude of the dipole moment. The factor cos2i5m, expresses the degree of po- 
larization: in a strong field, 6^ = and the dipole is at maximum strength, 
whereas in zero field Sm = 7r/4 and the dipole vanishes. Notice that for a given 
value of TO, the potential generated by the states with e = ± differ by a sign, 
this is appropriate, since these states correspond to dipoles pointing in opposite 
directions (Figure 1). 

The off-diagonal matrix elements in (PT|) are also important, for two reasons. 
First, it may be desirable to create superpositions of different energy eigenstates, 
and computing these matrix elements requires the off-diagonal elements in (I2ip . 
as we will see shortly. Second, when two dipoles interact with each other, one 
will experience the electric field due to the other, and this field need not lie 
parallel to the z-axis. Hence, the to quantum number of an individual dipole is 
no longer conserved, and elements of (j2ip with q =/= are required. 

4.1 Example: j = 1/2 

To illustrate these abstract points, we consider here the simplest molecular state 
with a A doublet: a molecule with j = 1/2, which has oj — 1/2, and consists of 
four internal states in our model. Based on the discussion above, we tabulate 
the matrix elements between these states in Table I. Because we have j = 1/2, 
we suppress the index j in this section. 

For concreteness, we focus on a type \a) molecule, as defined in Figure 1. 
This molecule aligns with the field and produces an electrostatic potential ([^5]) . 
However, if the molecule is prepared in a state that is a superposition of this 
state with another, a different electrostatic potential can result. We first note 
that combining \a) with \b) produces nothing new, since both states generate 
the same potential. 

An alternative superposition combines states \a) with state |c). In this case 
the two states have the same value of to, but are nevertheless non-degenerate. 
We define 

|V>),, = +) + Be— "^licS, -), 

for arbitrary complex numbers A and B with + = 1. Because the two 
states are non-degenerate, it is necessary to include the explicit time-dependent 
phase factors, where ujq = \Emii;e\/^, and 2hujo is the energy difference between 
the states. As usual in quantum mechanics, these phases will beat against one 
another to make the observables time dependent. 

Now some algebra identifies the mean value of the electrostatic potential, 
averaged over state |V')ac, as 

ac{i^mm)ac (24) 

= [(I - COS 2Sif2 -2\AB\ sin 2(5i/2 cos(2woi - S)] cos 0. 

This potential has the usual cos 9 angular dependence, meaning that the dipole 
remains aligned along the field's axis. However, the magnitude, and even the 
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Table 1: Matrix elements {mLue\Ciq\m'oje') for a j = 1/2 molecule. To obtain 
matrix elements of the electrostatic potential these matrix elements should 
be multiplied by (—l)'^Ck-q{d(j))fJ,/r'^, where q = m — m'. 









1 - ^'^+) 


1 - ^^-> 




i cos2(5i/2 


sin2(5i/2 


3 







— i sin 2(^1/9 


-icos2(5i/2 





3 




3 





-i COS 2^1/2 


-i sin2(5i/2 







3 


sin25i/2 


i cos2(5i/2 



sign, of the dipole change over time. The first term in square brackets in (|24p 
gives a constant, dc component to the dipole moment, which depends on the 
population imbalance |v4p — between the two states. The second term adds 
to this an oscillating component with angular frequency 2luo. The leftover phase 
5 is an irrelevant offset, and comes from the phase of A*B, i.e., the relative phase 
of the two components at time t = 0. 

It is therefore possible to construct a superposition of states of the dipole, 
such that the effective dipole moment of the molecule bobs up and down in time. 
The amount that the dipole bobs, relative to the constant component, can be 
controlled by the relative population in the two states. Moreover, the degree 
of polarization of the molecule plays a significant role. For a fully polarized 
molecule, when S1/2 = 0, only the dc portion of the dipole persists, although 
even it can vanish if there is equal population in the two states, "dipole up" and 
"dipole down." 

As another example, we consider the superposition of \a) with \d). Now the 
two states have different values of m as well as different energies: 



ad 



'2 ' ' '2 
The dipole potential this superposition generates is 

= JL cos 2^1/2 (I A|2 - \Bn COS e + I^Re U*Be'*(^+^'-''*) 

This expression can be put in a useful and interesting form if we parametrize 
the coefficients A and B as 

2 2 

This way of writing A and B seems arbitrary, but it is not. It is the same 
parametrization that is used in constructing the Bloch sphere, which is a pow- 
erful tool in the analysis of any two- level system [6]. 



16 



This parametrization leads to the following expression for the potential: 

ad{M'^m)ad (25) 

= --^ [cos2(5i/2 cos a cos 6* + sinQ;sin6'cos((/3 — 2LOot) — 0)] . 

The interpretation of this result is clear upon comparing it to the classical result 
First consider that the molecule is perfectly polarized, so that cos2(5i/2 — 1. 
Then ((25)) represents the potential due to a dipole whose polar coordinates are 
— ojt). That is, this dipole makes (on average) an angle a with respect 
to the field, and it precesses about the field with an angular frequency 2wo- 
Interestingly, even in this strong field limit where the field nominally aligns the 
dipole along z, quantum mechanics allows the dipole to point in quite a different 
direction. As the field relaxes, the z-component reduces, but this dipole still 
has a component precessing about the field. 



4.2 Example: j = 1 

We also consider a molecule with spin j — 1. Here there are in principle three 
mixing angles, (5i, (5o, and However, as noted above we have = —5i and 
(^o = T^l^i so that the entire electric field dependence of these matrix elements is 
incorporated in the single parameter 5i. In this notation, the matrix elements 
of the electrostatic potential for a, j — 1 molecule are given in Table 2. 

Similar remarks apply to the spin-1 case as applied to the spin-1/2 case. If 
the molecule is in an eigenstate, say | + Iw,— ), then the expectation value of 
the dipole points along the field axis, and its distribution has the usual cosO 
dependence. In an eigenstate with to = 0, however, the expectation value of the 
dipole vanishes altogether. 

As before, the molecule can also be in a superposition state. No matter 
how complicated this superposition is, the expectation value of the dipole must 
instantaneously point in some direction, since the only available angular depen- 
dence resides in the Ciq functions, which yield only dipoles. In other words, 
no superposition can generate the field pattern of a quadrupole moment, for 
example. 

Where this dipole points, and how its orientation evolves with time, however, 
can be complicated. For example, a superposition of | + Iw, +) and | + Iw, — ) can 
bob up and down, just like the analogous superposition for j = 1/2. However, 
for j = 1 molecules additional superpositions are possible. For example, consider 
the combination 

|V')3 = Ae"^°'\ + Iw, +> + Be''^^*|Ow, +) + Ce*"°*| - Ilu, -), 

where ?luia = A/2 is a shorthand notation for half the lambda doubling energy. 
Let us further assume for convenience that A, B, and C are all real. Then the 
expectation value of the electrostatic potential is 

3(V^|$(r)|^)3 = |cos25i(A2 + C^)2^ 
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+ ^cos(5i+7r/4)S^ 

X [Acos{{u)A — wo)f — <p) — Ccos(— (wa — i^o)t — (f>)] ■ 

By analogy with remarks in the previous section, this represents a dipole with 
a constant component along the 0-axis, which depends on both the strength 
of the field and on |Ap + |Cp, the total population in the ±to states. It also 
has a component in the x-y plane, orthogonal to the field's direction. In the 
case where C = 0, this component would precess around the field axis with a 
frequency to a — wq, in a clockwise direction as viewed from the +z direction. 
Vice-versa, if ^ = 0, this component would rotate at this frequency but in a 
counter-clockwise direction. If both components are present and A = C, then 
the result will be, not a rotation, but an oscillation of this component from, say, 
+x to —X, in much the same way that linearly polarized light in a superposition 
of left- and right-circularly polarized components. More generally, if A 7^ C, 
then the tip of the dipole moment will trace out an elliptical path. However, in 
the limit of zero field, uq reduces to oja and these time- dependent effects go 
away. 
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Table 2: Matrix elements of Ciq for a j = 1 molecule. To obtain the matrix elements of the electrostatic potential ^>(r), these 
matrix elements should be multiplied by {—l)'^Ck-q{04>)iJ./r^, where q = m — m'. 
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1 + Ico-) 
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{+lu)~\ 
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(-Iw+I 


^ cos 2^1 
— ^ sin 2di 
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-icos(5i +7r/4) 
i sin(<5i + 7r/4) 



1 cos(-(5i + 7r/4) 


i sin((5i + 7r/4) 
1 cos(5i + 7r/4) 





-isin(-5i +7r/4) 
-|cos(-5i +7r/4) 





-icos(-(5i +7r/4) 
i sin(-(5i + 7r/4) 

— 1 cos2(5i 

— 1 sin 2^1 





i sin(-5i + 7r/4) 
|cos(-(5i +7r/4) 
— ^ sin 26i 


{-lu!-\ 








-isin(-5i +7r/4) 


\ cos 2^1 



5 Interaction of dipoles 



Having thus carefully treated individual dipoles and their quantum mechanical 
matrix elements, we are now in a position to do the same for the dipole-dipole 
interaction between two molecules. This interaction depends on the orientation 
of each dipole, fli and fl2', and on their relative location, R. This interaction 
has the form (Appendix A) 



i?3 



i?3 



In going from the first line to the second, we assume that both molecules have 
the same size dipole moment fj,, and that the intermolecular axis makes an angle 
9 with respect to the laboratory z axis, so that R — {R, 0, (j)). The angles and 
(j) thus stand for something slightly different than in the previous section. The 
third line in (j26|) rewrites the interaction in a compact tensor notation that is 
useful for the calculations we are about to do. Here 

[Ml ® M2]2, = VS^I-I)" ( ^ \ Ci,,(/3iai)Ci,,(/?2a2) 

9l>32 ^ ^ 

denotes the second-rank tensor composed of the two first-rank tensors (i.e., 
vectors) Cigj^ (/3iQ!i) and Ciq^{f32Ci2) that give the orientation of the molecular 
axes [T]. Equation (pS)) highlights the important point that the orientations 
of the dipoles are intimately tied to the relative motion of the dipoles: if a 
molecule changes its internal state and sheds angular momentum, that angular 
momentum may appear in the orbital motion of the molecules around each 
other. 



5.1 Potential matrix elements 

Equation ((26|) is a perfectly reasonable way of writing the classical dipole-dipole 
interaction. Quantum mechanically, however, we are interested in molecules 
that are in particular quantum states \jmuj,t), rather than molecules whose 
dipoles point in particular directions [a, (3). We must therefore construct matrix 
elements of the interaction potential ([26| in the basis we have described in Sec. 
3.3. 

Writing the interaction in the form above has the advantage that each term in 
the sum factors into three pieces: one depending on the coordinates of molecule 
1, another depending on the coordinates of molecule 2, and a third depending on 
the relative coordinates (0, 0). This makes it easier to evaluate the Hamiltonian 
in a given basis. For two molecules we consider the basis functions 

(Q;i/3i|jmiw, ei){a2l32\jm2U), 62), (27) 
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as defined above. In this basis, matrix elements of the interaction become 

(jmiw, ei; jTO2W, £2! 14(6*, (j))\jm[id, e^^jm'^u}, e'2) = 
\/30/^2 / 2 1 1 



i?3 y q -q^ 
X (jTOiW, ei\Ciq, \jm[ui, e[) {jm2Ui, e2\Ciq2 Ijm^ui, e2)C2_g( 



(28) 



where matrix elements of the form (jmw, e|Cig|jm'Lj, e') are evaluated in Eq. 
(j22p . Conservation of angular momentum projection constrains the values of 
the summation indices, so that qi = mi — m'^, q2 — m2 — m^, and q — qi + q2 
= (mi + m2) — (m'l + m'2). To make this model concrete, we report here the 
values of the second-rank reduced spherical harmonics [1] : 

C20 = i(3cos2(?-l) 
/o\ 1/2 

C2±i = ^[2) cosesinOe^"!' 



1/2 



,2 /)„±2i<< 



C2±2 = y-^j sin- 

We also tabulate the relevant i-j symbols in Table III. 

Viewed roughly as a collision process, we can think of two molecules ap- 
proaching each other with angular momenta mi and m2, scattering, and de- 
parting with angular momenta m'l and m'2, in which case q is the angular mo- 
mentum transferred to the relative angular momentum of the pair of molecules. 
Remarkably, apart from a numerical factor that can be easily calculated, the 
part of the quantum mechanical dipole-dipole interaction corresponding to an- 
gular momentum transfer q has an angular dependence given simply by the 
multipole term C2-q- 

Suppose that the molecules, when far apart, are in the well-defined states 
((27|) . Then the diagonal matrix element of the dipole-dipole potential evaluates 
to 

miciili \ ( m2e2W ^ \ (l-3cos^' 



This has exactly the form of the interaction for classical, polarized dipoles, as 
in Eq. fS]). The difference is that each dipole /j, is replaced by a quantum- 
corrected version (in large parentheses). It is no coincidence that this is the 
same quantum-corrected dipole moment that appeared in the expression (|23p 
for the field due to a single dipole. When both dipoles are aligned with the 
field, we have miei > and m2e2 > (e.g., both molecules are of type |a)), 
and the interaction has the angular dependence oc (1 — 3cos^ 0). On the other 
hand, when one dipole is aligned with the field and the other is against (e.g., 
one molecule is of type \a) and the other is of type |c)), then the opposite sign 
occurs - just as we would expect from classical intuition. 

More generally, at finite electric field, or at finite values of R, the molecules do 
not remain in the separated- molecule eigenstates (|27p , since they exert torques 
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Table 3: The 3-j symbols needed to construct the matrix elements in (^5]) . 
Note that these symbols remain invariant under interchanging the indices qi 
and (72, as well as under simultaneously changing the signs of g — 1, q2, and q 
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-1 


l/\/30 
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1/V5 



on one another. The interaction among several different internal molecular 
states makes the scattering of two molecules a "multichannel problem," the 
formulation and solution of which is described in Chapters XXX. However, a 
good way to visualize the action of the dipole-dipole potential on the molecules 
is to construct an adiabatic surface. To do so, we diagonalize the interaction at 
a fixed value of the relative location of the two molecules. 

Before doing this, we must consider the quantum statistics of the molecules. 
If the two molecules under consideration are identical bosons or identical fermions, 
then the total two-molecule wave function must account for this fact. This total 
wave function is 

(alA7l|J"^l^I^el)(a2/3272)b'n^-2'^e2)-P,■^^;mlelm2e2(^:^)<^)• 

This wave function is either symmetric or antisymmetric under the exchange of 
the two particles, which is accomplished by swapping the internal states of the 
molecules, while simultaneously exchanging their center-or-mass coordinates, 
i.e., by mapping Rio —R: 

(q;i/3i) ^ {oL2(ii) 

R ^ R (30) 

e TT-9 

(f) — > TT + (f). 

For the molecule's internal coordinates, a wave function with definite exchange 
symmetry is given by 

(q;i/3i7i I jmiwei) (02/3272) | j'7i2<De2)s 



■\/2(l + 6mim2^eie2) 

X [{aiPiji\jmiUjei){a2f32j2\jm2Uje2) + s{a2P2l2\jmiUjei) {aiPiji\jm2Uje2)] ■ 

The new index s = ±1 denotes whether the combination (|3ip is even or odd 
under the interchange. If s = +1, then F must be symmetric under the trans- 
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formation R — > —R for bosons, and odd under this transformation for identical 
fermions. If s — —1, the reverse must hold. 

We now have the tools required to consider the form of the dipole-dipole 
interaction beyond the "pure" dipolar form P5)) . The details of this analysis 
will depend on the Schrodinger equation to be solved. In its fundamental form, 
the Schrodinger equation reads 



2m, 



Here Hs stands for the threshold Hamiltonian that includes A-doubling and 
electric field interactions, and is assumed to be diagonal in the basis (I5T|) : and 
rrir is the reduced mass of the pair of molecules. In the usual way, we expand 
the total wave function ip as 



where the index i stands for the collective set of quantum numbers {juj', mieim2t2s\ ■ 

Inserting this expansion into the Schrodinger equation and projecting onto 
the ket {i\ leads to the following set of coupled equations: 



E 



dR^ 



1 



d 



R^ sin 6 80 X'"''' 80 



8 



1 



92 



i?2 sin^ ' 



Fi 



If we keep N channels i, then this represents a set of N coupled differential equa- 
tions. We can, in principle, solve these subject to physical boundary conditions 
for any bound or scattering problem at hand. For visualization, however, we 
will find it convenient to reduce these equations to fewer than three independent 
variables. We carry out this task in the following subsections. 



5.2 Adiabatic potential energy surfaces in two dimensions 

Applying an electric field in the z direction establishes z as an axis of cylindrical 
symmetry for the two-body interaction. The angle determines the relative 
orientation of the two molecules about this axis, thus the interaction cannot 
depend on this angle. To handle this, we include an additional factor in our 
basis set, 

\mi) = -^=cyip{imi(f>). 
V 27r 

We then expand the total wave function as 
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In each term of this expression the quantum numbers must satisfy the con- 
servation requirement for fixed total angular momentum projection, Mtot = 
mi + 7712 +mi. In addition, applying exchange symmetry to each term requires 
that (i?, TT - 61) = s(-l)"'i^,,™,(i?,6') for bosons, and s(-l)"'+iF,,„, (i?, 6*) 
for fermions. 

Inserting this expansion into the Schrodinger equation yields a slightly dif- 
ferent set of coupled equations: 



2mr 



1 d f . „d 
sini 



Ti m. 



2mrR'^ sin^ 9 



F^rn,+J2('\Vi''\i')F,,„,, + (^\H s\i) F,„,^ = EF, 



mil ■ 



This substitution has the effect of replacing the differential form of the azimuthal 
kinetic energy, c>c d'^/d4>^, by an effective centrifugal potential (x m^/R^svo? 6. 
In addition, the matrix elements of the dipolar potential V]l^ are slightly dif- 
ferent from those of Vd- Recall that the 9 dependent part of the matrix element 
(|28p is proportional to C2-q{9, <j)), which we will write explicitly as 

C2^q{9,^) = C2-q{9)cM-iq(t>). 

This equation explicitly defines a new function C2-q{9) that is a function of 
6 alone, and that is proportional to an associated Legendre polynomial. The 
matrix element of the potential now includes the following integral: 

{mi\C2-q\m[) = ( d0-^e-™"^C2-9(^)e-"?'^-^e™2'^ 
J v27r \/27r 



= ^M,^,M'^^C2-q{9), 



which establishes the conservation of the projection of total angular momentum 
by the dipole-dipole interaction. Therefore, matrix elements of V^^ in this 
representation are identical to those in of Vd in Eq. (|28p except that the factor 
exp(-ig(/)) is replaced by Sm^^^m^^^- 

With these matrix elements in hand, we can construct solutions to the cou- 
pled differential equations ([5^ . However, to understand the character of the 
potential surface, it is useful to construct adiabatic potential energy surfaces. 
This means that, for a fixed relative position of the molecules {R,9), we find 
the energy spectrum of ([5^ by diagonalizing the Hamiltonian V^^ + V^^ + Hs, 
where V^^ is a shorthand notation for the centrifugal potential discussed above. 
This approximation is common throughout atomic and molecular physics, and 
amounts to defining a single surface that comes as close as possible to repre- 
senting what is, ultimately, multichannel dynamics. 
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5.3 Example: j = 1/2 molecules 



Analytic results for the adiabatic surfaces are rather difficult to obtain. Consider 
the simplest realization of our model, a molecule with spin j — 1/2. In this 
case each molecule has four internal states (two values of m and two values of 
e), so that the two-molecule basis comprises sixteen elements. Dividing these 
according to exchange symmetry of the molecules' internal coordinates, there 
are ten channels within the manifold of s = +1 channels, and six within the 
s — ~1 manifold. These are the cases we will discuss in the following, although 
the same qualitative features also appear in higher-j molecules. 

As the simplest illustration of the influence of internal structure on the dipo- 
lar interaction, we will focus on the lowest-energy adiabatic surface, and show 
how it differs from the "pure dipolar" result l|29|) as the molecules approach one 
another. The physics underlying this difference arises from the fact that the 
dipole-dipole interaction becomes stronger as the molecules get closer together, 
and at some point this interaction is stronger than the action of the external 
field that holds their orientation fixed in the lab. The intermolecular distance 
at which this happens can be approximately calculated by setting the two in- 
teractions equal, /^^/i?o ~ V (/^^)^ + yielding a characteristic distance 



When Rq, the electric field interaction is dominant, the dipoles are aligned, 
and the interaction is given by Eq. ((29|) . When R becomes comparable to, or 
less than, i?Oj then the dipoles tend to align in a head-to-tail orientation to 
minimize their energy, regardless of their relative location. 

Before proceeding, it is instructive to point out how large the scale Rq can 
be for realistic molecules. For the OH molecule considered above, with = 1.7 
Debye and A = 0.06 cm~^, the molecule can be polarized in a field of £ w 1600 
V/cm. At this field, the characteristic radius is approximately Rq w 120 oq 
(where qq = 0.053 nm is the Bohr radius), far larger than the scale of the 
molecules themselves. Therefore, while the dipole-dipole interaction is by far 
the largest interaction energy at large R, over a significant range of this potential 
does not take the usual dipolar form. To take an even more extreme case, the 
molecule NiH has a ground state of symmetry with j = 5/2 [2]. Because 
it is a A, rather than a H, state, its A doublet is far smaller, probably on the 
order of ~ 10~^ cm~^. This translates into a critical field of £ w 0.5 V/cm, 
and a characteristic radius at this field Rq « 2000 gq « 0.09 //m. This length 
is approaching a non-negligible fraction of the interparticle distance in a Bose- 
Einstein condensed sample of such molecules (assuming a density of 10^^ cm~^, 
this spacing is of order 0.2 /im). Deviations from the simple dipolar behavior 
may thus influence the macroscopic properties of a quantum degenerate dipolar 
gas. 

As an example, we present in Figure 2 sections of the lowest-energy adiabatic 
potential energy surfaces for a ficticious j ~ 1/2 molecule whose mass, dipole 
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moment, and A doublet are equal to those of OH. These were calculated in a 
strong-field limit with £ = 10^ V/cm. Each row corresponds to a particular 
intermolecular separation, which is compared to the characteristic radius Rq. 
However, as noted above, there are two possible ways for the molecule to have 
its lowest energy, as illustrated by parts (a) and (b) of Figure 1. Interestingly, it 
turns out that these give rise to rather different adiabatic surfaces. To illustrate 
this, we show in the left column of Figure 2 the surface for a pair of type \a) 
molecules, which corresponds at infinitely large R to the channel \\ + , \+] s = 
1); and in the right column we show combinations of one type \a) and one type 
1 6) molecule. In the latter case, there are two possible symmetries corresponding 
to s = ±1, both of which are shown. Finally, for comparison, the unperturbed 
"pure dipole" result is shown in all panels as a dotted line. 

Consider first two molecules of type |a) (left column of Fig. 2). For large 
distances R > Rq (top panel), the adiabatic potential deviates only slightly from 
the pure-state result, reducing the repulsion at = 7r/2. When R approaches 
the characteristic radius Rq (middle panel), the effect of mixing in the higher- 
energy channels becomes apparent. Finally, when R < Rq (lower panel), the 
mixing is even more significant. In this case the dipole-dipole interaction is the 
dominant energy, with the threshold energies serving as a small perturbation. 
As a consequence, the quantum numbers s = 1) can no longer identify 

the channel. It is beyond the scope of this chapter to discuss the corresponding 
eigenstates in detail. Nevertheless, we find that for R < Rq the channel \aa) 
(repulsive at = 7r/2) is strongly mixed with the channel \ac) (attractive at 
9 = tt/2). The combination is just sufficient that the two channels nearly 
cancel out one another's ^-dependence. At the ends of the range, however, 
the adiabatic curve is contaminated by a small amount of channels containing 
centrifugal energy oc 1/ sin^ 9. 

The right column of Figure 2 shows adiabatic curves for the mixed channels, 
one molecule of type \a) and one of type \b). In this case there are two possible 
signs of s; These are distinguished by using solid lines for channel \^+,~^~; s = 
1); and dashed lines for channel |^+,— ^ — ;s — —1). Strikingly, these surfaces 
are different both from one another, and from the surfaces in the left column 
of the figure. Ultimately this arises from different kinds of channel couplings 
in the potentials l|28|) . Note that, while type \a) and type \b) molecules are 
identical in their interaction energy with the electric field, they still represent 
different angular momentum states. Nevertheless, molecules in these channels 
still closely reflect the pure dipolar potential at large R, and become nearly 
^-independent for small R. 

A further important point is that the potentials described here represent 
large energies as compared to the mK or fiK translational kinetic energies of cold 
molecules, and will therefore significantly influence their dynamics. Further, the 
potentials depend strongly on the value of the electric fleld of the environment, 
both through the direct effect of polarization on the magnitude of the dipole 
moments, and through the influence of the flcld on the characteristic radius Rq. 
It is this sensitivity to field that opens the possibility of control over interactions 
in an ultracold dipolar gas. 
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Figure 2: Angular dependence of adiabatic potential energies for various com- 
binations of molecules at different intcrparticle spacings i?,, which arc indicted 
on the right side. Dotted lines: diagonal matrix element of the interaction, 
assuming both molecules remain strongly aligned with the electric field. Solid 
and dashed lines: adiabatic surfaces. These surfaces are based on the j = 1/2 
model discussed in the text, using fi = 1.68 Debye, A = 0.056 cm~^, = 8.5 
amu, and £ = 10* V/cm, yielding Rq ~ 70ao- The left hand column presents 
results for two molecules of type |a), as labeled in Fig. 1; in the right column 
are results for one molecule of type \a) and one of type \b), which necessitates 
specifying an exchange symmetry s. 

0.02: J^^^ , 



^ 0.00 



-0.02 
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Although we have Hmited the discussion here to the lowest adiabatic state, 
interesting phenomena are also expected to arise due to avoided crossings in 
excited states. Notable is a collection of long-range quasi-bound states, whose 
intermolecular spacing is roughly centered around Rq [7] . Such states could con- 
ceivably be used to associate pairs of molecules into well-characterized transient 
states, furthering the possibilities of control of molecular interactions. 

5.4 Adiabatic potential energy curves in one dimension: 
partial waves 

For many scattering applications, it is not necessarily convenient to express 
the dipole-dipole interaction as a surface (more properly, a set of surfaces) in 
the variables {R, 9, (p) describing the relative position and orientation of the 
molecules. Rather, it is useful to expand the relative angular coordinates {6, (j)) 
in a basis as well. To do this, the basis set ([?7|) is augmented by spherical 
harmonics describing the relative orientation of the molecules, to become 

(ai/3i7i I jmiwei) (a2^272 \jm2UJe2){0(i)\lmi) , 

with 



{e<j)\lmi) = Yi„,,{B<j)) = ^——Ci,n^{B<j)). 
The total wave function is therefore described by the superposition 

i' .V .m'^ 

which represent a conventional expansion into partial waves. The wave function 
is, as above, restricted by conservation of angular momentum to require Mtot = 
mi + mi + 771/ to have a constant value. Moreover, the effect of the symmetry 
operation (jSO]) on \lmi) is to introduce a phase factor (—1)'. Therefore, the wave 
function is restricted to s(— 1)' = 1 for bosons, and s(— 1)' = —1 for fermions. 

The effect of this extra basis function is to replace the C2-q factor in (|28p 
by its matrix element 

C2_, ^ {lmi\C2^q\l'm[) (33) 

= V(2^ + 1)(2.-M)(-1)-(J ™0 

From this expression it is seen that the angular momentum g, lost to the internal 
degrees of freedom of the molecules, appears as the change in their relative 
orbital angular momentum. From here, the effects of field dressing are exactly 
as treated above. 

The quantum number Z, as is usual in quantum mechanics when treated in 
spherical coordinates, represents the orbital angular momentum of the pair of 
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molecules about their center of mass. Following the usual treatment, this leads 
to a set of coupled radial Schrodinger equations for the relative motion of the 
molecules: 

2mr dR^ 2mrR^ 

i' 

The second term in this expression represents a centrifugal potential cx 
which is present for all partial waves ^ > 0. 



5.5 Asymptotic form of the interaction 

Casting the Schrodinger equation as an expansion in partial waves, and the 
interaction as a set of curves in R, rather than as a surface in (i?, 0, </)), allows us 
to explore more readily the long-range behavior of the dipole-dipole interaction. 

The first 3-j symbol in p3)) vanishes unless ^ + 2 + T is an odd number, 
meaning that even(odd) partial waves are coupled only to even(odd) partial 
waves. Moreover, the values of I and I' can differ by at most two. Thus the 
dipolar interaction can change the orbital angular momentum state from I = 2 
to /' = 4, for instance, but not to V — 6. Finally, the interaction vanishes 
altogether for I = V = Q, meaning that the dipolc-dipolc interaction nominally 
vanishes in the s-wave channel. 

Since all other channels have higher energy, due to their centrifugal poten- 
tials, it appears that the lowest adiabatic curve is trivially equal to zero. This 
is not the case, however, since the s-wave channel is coupled to a nearby d- 
wave channel with I = 2. Ignoring for the moment higher partial waves, the 
Hamiltonian corresponding to a particular channel i at long range has the form 

/ Ao2/i?3 \ 

A2o/i?^ -if? IrUrR^ ^ A22IR? ) ■ 

Here Aqi = ^20 and ^22 are coupling coefficients that follow from the ex- 
pressions derived above. Note that all these coefficients are functions of the 
electric field E. Now, the comparison of dipolar and centrifugal energies de- 
fines another typical length scale for the interaction, namely, the one where 
[1^ jR? = /rrirR^, defining a "dipole radius" Rd — ^J^rrir/h^ . (More properly, 
one could define an electric-field-dependent radius by substituting ^02 for /i^). 
For OH, this length is 6800 ao, while for NiH it is ^ 9000 ao. 

When the molecules are far apart, R > Rd, the dipolar interaction is a 
perturbation. The size of this perturbation on the s-wave interaction is found 
through second-order perturbation theory to be 

m^/mrR^ ^ \ ) R^' 
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Therefore, at very large intermolecular distances, the effective potential, as de- 
scribed by the lowest adiabatic curve, carries a dependence on R, rather 
than the nominal dependence. At closer range, when R < Rjj, the 1/R^ 
terms dominate the 1 /i?^ centrifugal interaction, and the potential reduces again 
to the expected 1/R^ dependence on R. 
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